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Abstract. Recent experimental results from supernovae la observations have been 
interpreted to show that the rate of expansion of the universe is increasing. Other 
recent experimental results find strong indications that the universe is “flat.” In this 
paper, I investigate some solutions of Einstein’s field equations which go smoothly 
between Schwarzschild’s relativistic gravitational solution near a mass concentration 
to the Friedmann-Lemaitre expanding universe solution. In particular, the static, 
curved-space extension of the Lemaitre-Schwarzschild solution in vacuum is given. 
Uniqueness conditions are discussed. One of these metrics preserves the “cosmologi¬ 
cal equation.” We find that when the rate of expansion of the universe is increasing, 
space is broken up into domains of attraction. Outside a domain of attraction, the 
expansion of the universe is strong enough to accelerate a test particle away from the 
domain boundary. I give a domain-size-mass relationship. This relationship may very 
well be important to our understanding of the large scale structure of the universe. 

pacs: 96.30Sf, 97.60Bw, 95.85Bh, 04.20-q 


1. INTRODUCTION AND SUMMARY 


Recently de Bernardis et al.[ 1] reported that a very careful examination 
of the cosmic microwave background provides, as further explained by 
Hu[2], very strong evidence that the universe is flat! 

In addition, data has recently been reported by Riess et al.[ 3], and 
Perlmutter et al.[ 4] on measurements of the luminosity and the redshift 
of a number of high-redshift supernovae of type la. The authors’ best 
fit to their data involves a strongly curved space-time. The curvature 
constant (defined below) that they find is about -1.05, which corre¬ 
sponds to a strongly curved, open universe. However, their error bars 
are sufficiently large so that a flat universe is not inconsistent with their 
data. These authors have analyzed their data on the basis of the “cos¬ 
mological equation” between the rate of expansion, the mean energy 
density, the radius of curvature of space, and the cosmological con¬ 
stant. This relation has been derived from Einstein’s field equations[5] 
under the assumption of the validity of the Friedmann-Lemaitre line 
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element at large distances. The authors have concluded, in terms of this 
model of the universe, that instead of the rate of expansion decreasing, 
as many workers had thought, their data is best fit by a model in 
which the rate of expansion is increasing. Riess et al.[ 3] find that the 
deceleration/acceleration parameter (as defined in section VI) qo < 0 
(acceleration) with better than 90% confidence, and Perlmutter et al. [4] 
find the same at the 2.86 standard deviation level of confidence. 

In this paper, I investigate some different metrics which may fa¬ 
cilitate the investigation of some of the consequences of the reported 
acceleration of the expansion of the universe. I will focus on spherically 
symmetric models with a central mass concentration. 

The fact that the rate of expansion of the universe is increasing 
leads to the conclusion that the universe is broken up into domains of 
attraction. Briefly, the underlying physics of this feature may be seen in 
the following rough analysis. The equation of purely radial motion of a 
test particle at rest with respect to the Friedmann-Lemaitre coordinate 
system is just 

r = -r, (1) 


where r is the proper distance, and a is the universal expansion factor 
which is a function of time alone. The most simple way to approximate 
the effects of a gravitating mass concentration is just to add Newton’s 
force term so that 


a 

f = —r — 
a 


GM 
r 2 ' 


( 2 ) 


If the expansion of the universe is decelerating, then a < 0, so r < 0 
always. Near a gravitating mass the expansion is clearly unimportant 
and Newton’s laws hold with only minimal corrections. On the other 
hand, if the expansion is accelerating, then a > 0. Thus there will be 
a distance such that the expansion of the universe exactly balances 
the gravitational attraction. Test particles at smaller distances will be 
accelerated inwards. I call this region a domain of attraction. Test parti¬ 
cles outside will be accelerated outwards. Hence in this very simple case 
there are two attractors, one is the mass concentration and the other is 
the point at infinity. We give a domain-size-mass relation in equations 
(73) and (76). In the decelerating case, there is just one attractor and 
its domain of attraction is the entire space. 

In the second section we remind the reader of the Lemaitre-Tolman 
formalism. This formalism has the feature that the proper time appears 
explicitly. Thus space-time is divided into a set of three-dimensional, 
space-like manifolds which are indexed by the proper time. Since it is 
the proper time which appears in the universal expansion function in 
the Friedmann-Lemaitre line element and it governs the universe as a 
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whole, it seems important to preserve the property that the universal 
expansion factor is a function of the proper time alone. 

In the third section we review three different ways to go from the 
Schwarzschild metric of a planetary system to expanding space at large 
distances. The McVittie solution has a problem at the Schwarzschild ra¬ 
dius using normal expansion factors. The Einstein-Straus “Swiss Cheese 
Model” is unstable to perturbations and there are orbits which are dis¬ 
continuous functions of their initial conditions. The Bona-Stela model 
using Liebovitz metric insertions in a Friedmann-Lemaitre background 
predicts an increase in the length of the earth’s year which is in strong 
disagreement with observations. 

In the fourth section we consider the idea that the vacuum is some¬ 
thing more than just empty space. The self-energy of the vacuum 
may correspond to a mass-energy density of empty space. We know 
on the laboratory scale from, for example, the Casimir effect, and 
on the microscopic scale from particle physics that the vacuum does 
have measurable effects. So far as I know, particle theorists have not 
yet been able to compute a quantitative result for the energy density. 
However, they do feel that such an idea is extremely plausible. We 
consider in this section a solution to Einstein’s held equations which 
has a homogeneous vacuum, mass-energy density and so preserves the 
so called “cosmological equation.” This solution reduces precisely to 
the exterior Schwarzschild solution near the central mass condensa¬ 
tion. In addition, it reduces to the Friedmann-Lemaitre solution far 
from the mass concentration. This solution is continuously, infinitely 
differentiable everywhere, except at the mass concentration. 

In the fifth section, I explore some of the properties of the solution 
obtained in the fourth section. Both “hat” and curved space are consid¬ 
ered. The limiting results are as expected and the transition between 
the two aforementioned limits are illustrated. In particular the static, 
curved-space extension of the Lemaitre-Schwarzschild solution is given 
for a mass concentration in a vacuum. 

In the sixth section, I consider the implications of the increasing 
rate of expansion of the universe on some of the large scale structures 
found in the universe. It is found that this feature creates domains 
of attraction. A domain-size-mass relation is derived. Outside these 
domains, the increasing rate of the expansion of the universe would, 
in time, be expected to tear the structures apart. It is suggested, from 
the correspondence between the predicted size of these domains of at¬ 
traction, the size of the Local Group, and the size of the Virgo Cluster, 
that this effect may well be important in any study of the large scale 
structure of the universe. 
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2. THE LEMAiTRE-TOLMAN FORMALISM 

I will find it convenient to employ the Lemaitre-Tolman[6, 7] formalism. 
In this case we start with the line element, 

ds 2 = -e 2a ^ r Up 2 - e 2 ^ p ' T \d6 2 + sin 2 0 d(j) 2 ) + c 2 dr 2 , (3) 

where c is the velocity of light, and r is the proper time for an ob¬ 
server at rest in this coordinate system. This line element is spherically 
symmetric in its structure, which is of sufficient generality for my 
purposes. 

The reason for our interest in this particular form of the line ele¬ 
ment is that two important line elements are special cases. First, the 
Friedmann-Lemaitre line element, 

ds 2 = — a 2 (r)dp 2 — a 2 (r)p 2 ( d6 2 + sin 2 6 dtp 2 ) + c 2 dr 2 (4) 

is manifestly of form (3). The second important line element is the 
static Schwarzschild line element. It is usually written in the form, 

ds 2 =-~~ ? ' 2 2 s ™ 2 ^ c ^ >2 ) + c 2 ^1 — df 2 (5) 

c 2 r 

However, as is well known[8] this form does not extend inside the 
Schwarzschild radius rs = 2GM/c 2 . An alternate form has been given 
by Lemaitre [9]. It is 

ds 2 = — ^J^ dp 2 — 7 Z 2 (d0 2 + sin 2 6 dtp 2 ) + c 2 dr 2 , 
r -| 2/3 

77 = %y/2GM/*(p-cr) . (6) 

This form, without going into the full generality of the Schwarzschild 
solution, has no singularity at the Schwarzschild radius. It is a coor¬ 
dinate system adapted to a freely falling observer. This metric can 
easily be seen to have the form (3). As it is our goal to find a metric 
which tends asymptotically to the Friedmann-Lemaitre metric at very 
large distances, and tends asymptotically to the Schwarzschild metric 
on the scale of our planetary system, it seems appropriate to use the 
Lemaitre-Tolman formalism. In this formalism, four dimensional space- 
time is described by a set of three-dimensional, space-like hyper surfaces, 
indexed by the time like variable r. 

Next we consider the stress-energy tensor for the line element (3). 
First, 

8ttcT f = -2 (V/3 - d/3' + /3') = 0, (7) 
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where an overdot means differentiation with respect to r and ' means 
differentiation with respect to p. The conditions Tf = —e 2 a T 4 = 0 are 
set in order to have a time-orthogonal coordinate system. The solution 
of this equation is well known to be, 


e 


a 


(d'e? 

mr 


0 < f(p) < 00 , 


( 8 ) 


where f(p) is undetermined. This result allows us to eliminate a(p,r) 
in terms of / 3(p,r) and f{p). 

The other non-zero elements of the stress-energy tensor are, 


8 ttTI 


8 tT| 


8vtT| 


e 2/3 — e 2 "(/3 ') 2 + 2 4 + 3 C-l “ A 


2 4 + 3 (-| + e~ 2/3 [l — f 2 (p)] — A 


(9) 


gwT , = j3 + P + a$ + a + a 2 _ ^ + {fj , ]2 _ ^ _ A 

+ C <>«) 
(7) + + e ~ 2p - e ~ 2a W+- ‘Zu'p'] - A 


p 

c 


+ 2 ^ + e~ 2p 
+ ^c 2 /3' +e 


1 - f\p) ~ 


mp)f(p) 


p' 


-A ( 11 ) 


We see that by (10), the necessary and sufficient condition that the 
spatial curvature be isotropic is that T ) 1 be independent of p. In addi¬ 
tion, this condition implies that the spatial curvature is homogeneous 
throughout all space. 

If we replace f(p) by unity, then by construction, e 13 = f e a dp, the 
distance measured from the origin. Hence the area of a sphere is just 47 r 
times the square of the radius, as is given by the standard Euclidean 
formula. I will call this special case the case of “flat” space. In the 
special case of “flat” space and isotropic spatial curvature, (9) and (11) 
become 


8 v tT} = 87tT 2 2 = 8 vrT 3 3 



( 12 ) 


8 ttT 4 4 


P 

c 




-A. 


(13) 
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We observe, that (9) and (11) may be rewritten as, 


d 

8t 


0 3/3 


P 

c 


A 

3 


+ e 0 [l- f 2 (p)} i = 8nT}Pe 30 


hr 



A 

3 


+ e 0 [l-f 2 (p)} 


8irT£/3'e 30 . 


The quantity 





(14) 

(15) 


(16) 


is easily recognized[6, 10] to be the mass equivalent to the total energy 
contained within the comoving shell with radial coordinate p at time 

T. 

We will also be concerned with the dynamics as seen in this for¬ 
malism. The main interest will be to assess the difference between the 
effects on planetary systems of the metrics which we will study herein, 
and the Schwarzschild metric. In particular I am interested in the solar 
system. The standard formula for the equations of motion of a test 
particle [5] is 


fds\ d 

\dr J dr 



1 dx l dxi 

2 9ij ’ k inr 


(17) 


where of course dx 4 /dr = 1. The fourth component of (17) gives the 
convenient equation, 




d-s 

dr 


-i 


ae 2a p 2 + /3e 20 \ 0 2 + sin 2 Q0k2 


Explicitly, the three equations of motion are, 


Be 2a p + — \e 2a p] = a'e 2a p 2 + (d'e 20 \o 2 + sin 2 9 ft 
dr L 

Be 20 O + — e 20 O = e 20 sin 9 cos 6 <j ) 2 , 
dr L J 


f d-s\ d 

\dr J dr 


e 20 sin 2 9 


ds 

dr 


-i 


= 0. 


(18) 

(19) 

( 20 ) 

( 21 ) 
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Following the completely standard modus operandi, we make an imme¬ 
diate simplification by setting 6 = n/2. Now equation (20) is explicitly 
satisfied and the other two equation are simplified. The first integral of 
(21) gives us the conservation of angular momentum, to wit, 


e 20 4> = C Q ( , 

\dr J 


( 22 ) 


where the Ci are constants of integration. By means of (3) we may 
write (22) as, 


<j> 2 = Cle-W (c 2 - e 2 V - e 2 ^ 2 ) = ^ 


de~ 4fj (c 2 - e 2 V) 


+ Cfe-W 


(23) 


These steps leave us with the single equation (19) for p(r) to deal 
with. It now becomes 


„ , , , 4> 2 Cfe-W ( c 2 - e 2 V) 

Bp+(2a + ap)p + p--- 

where B (18) may be re-expressed as, 

'p\ 2 CiPe-^ (c 2 - e 2 V) 


B = —ae 


2a 


c 2 (l + Cfie 2 ^) 


(24) 


(25) 


Thus, (19) becomes, 


p = -(2a+a'p)p-{- 


de- 2 ? (c 2 e~^ + PP'p) (c 2 - e 2a p 2 ) 
c 2 P' (1 + d e ~ W ) 


+pae 2a ' J 


(26) 

which gives an explicit, ordinary, non-linear differential equation for 
p(r). Coupled with 6 = vr/2 and equation (23), we have the equations 
of motion of the test particle. 

A more transparent form follows for the purpose of identifying the 
corrections to the Schwarzschild dynamics results if we make the change 
of variables, 

r = exp[p(p,r)], =*■ p = p(r,T). (27) 


Then we have 

Equation (19) becomes, 

*2 , /' (* 


f = [jd + P'pj r. 


r = 


p+p z 


fP’ \r 


~P 


c 2 d , p'+pp' 

r 3 \ c 2 f 2 P’ 


(28) 


■ \ 3 
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r 2 

L o 


+- 


— (r — fir) 
r 


r — /3r\ $(r — $r) 3 CqC 2 

rf J rc 2 f 2 r 4 

r (l + Cl/r 2 ) 


(29) 


3. CLASSICAL MODELS 


Currently there are two general relativistic descriptions of spacetime in 
popular use. For planetary systems and other gravitationally bound 
structures which are small on the scale of the universe, there is a 
static description of the behavior of spacetime. On the other hand, 
for large-scale behavior, there is a time dependent description which 
is appropriate as a description of phenomena such as the observed 
red-shift of distant galaxies. 

The classic question is, “How can these two disparate descriptions of 
spacetime possibly be reconciled with each other?” The first answer to 
this problem was given by McVittie[ll]. Using the line element form, 

ds 2 = — e M (dp 2 + p 2 dd 2 + p 2 sin 2 9d(j) 2 ) + e u dt 2 (30) 


He proposed the solution, 


= 


a(tf 


[1 + (a(t)p/2a(t)R) 


212 


1 + 


m \ 
2 a(t)p) 


e = c 


/1 — m/2a(t)p\ 2 

GM 
m = ——- 

1 + f—) 2 

\ 1 + m/2a(t)p ) 

c 2 

\2RJ 


1/2 


(31) 


where G is Newton’s constant of gravitation. Some of the virtues of 
this model can be seen by examining the corresponding stress-energy 
tensor, 


8v xTl 


8 ttT| 
8 t rT| 


1 + (m/2a(t.)p) 2 

a(t) 2 R 2 ( 1 + m/2ap) 5 (l — m/2a(t)p) 
2aa[l + m/2a(t)p] + a 2 [l — 5m/2a(t)p] 
a(t) 2 c 2 ( 1 — m/2a{t)p) 

8 ttT 3 3 = 8v tTI 

_3_ 3d 2 

a(t.) 2 R 2 (l + m/2a(t)p) 5 c 2 a 2 


In the limit as m —> 0 or p —> oo these tensor elements reduce to those 
corresponding to (4), and also in the limit as R —► oo and d —> 0 they 
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vanish as with the Schwarzschild metric. For this metric the spatial 
curvature is isotropic, but not homogeneous. The dynamical behavior 
of a test particle can be computed from (17). In the usual case 9 = n/2, 
the equation of motion, analogous to (29), becomes, 


r — 2 —r + r 
a 


1 , 


k-' 

2 -1 


1 a 


+ u r-r 

a 


JL fll LI -Li / I _ 

= + o (m ~ + 2 


a a* 


1 ... I 


A 2 a 4 e 


ds V 
dt) 


1 (ds\ z d 


—v ae ^ - 

2 2 


dt 7 dt 


-2 


a 

a 


(33) 


where, in concert with (22), r 2 e^(j) = A(ds/dt), and the change of 
variables r = a(t)p has been used. The slow speed limit of this equation 
of motion is, 

a A 2 c 2 GM 

r = -r H-;t-75- 

a r° 

This equation differs from Newton’s by the a term. This term is propor¬ 
tion to Hq where Ho ss 1.62 x 10” 18 /r,5o is Hubble’s constant. /150 = 1 if 
Hq = 50 knr/sec/Megaparsec. Thus in McVittie’s model, the effects on 
the solar system are too small to be currently measured. For systems 
like the solar system, the McVittie model works very well. However for 
other systems, if we make the popular assumption that a(t') oc t 7 then, 
as in McVittie’s solution a(t) is a function of coordinate time and not of 
proper time, there will be a singularity at the the Schwarzschild radius 
as a(t) will diverge there since t —► 00. 

A second answer to our question was given by Einstein and Straus 
[12]. It is the so called “Swiss Cheese Model.” They showed that the 
static Schwarzschild metric (5) can be matched to any Friednrann- 
Lemaitre-Robertson-Walker metric on a spherical surface of an expand¬ 
ing radius in the Schwarzschild metric, but of a constant radius in 
the Friedmann-Lemaitre-Robertson-Walker metric. It has been found 
however[13, 14] to be unstable with respect to perturbations. In ad¬ 
dition, there is another unphysical aspect of this model. The slow- 
speed dynamical equations in the interior Schwarzschild (r)) and in the 
exterior Friedmann-Lemaitre (r e ) regions are, respectively, 

d 2 r; GM _ d 2 f e d _ 

~dt 2 "' ~ 7 f r " ~d^ = a re 

For ease of exposition, we choose a{t ) = (t/to) 2 ^- Thus the general 
solution in the exterior region is r e = AG + Rt 1-7 = At 2 / 3 + Bt 1//3 , 
which is the equation of a parabola. 
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As an illustration of the behavior of the “Swiss cheese model” I have 
computed the following trajectories. I use as a unit of time the Hubble 
time, that is 1/Hq. As a unit of distance I use y/ MqG / Hq which is 
_ 9/3 

about 25 /i 50 ' million Astronomical units. Mq is the mass of the Sun. 
I set to = 1. The metric interface is at r = t 2 / 3 in these units for 
flat spacetime, as mentioned above. In order to follow a Friedmann- 
Lemaitre trajectory the test particle’s distance from the origin must be 
larger at every time than that for the interface. Let us take in rectan¬ 
gular coordinates the initial conditions, x = A, x = 0, y = 0, y = A. 
The trajectory in the external region is 



Thus, 

A 2 1 Of 4/3 - 22 1 + 13t 2/3 > f 4/3 . (37) 

The trajectory will intersect the interface if equation (37) is an equality. 
By means of the quadratic formula, an intersection will occur if 


t l/3 


11 ±^-9 + 13/A 2 
10-1/A 2 ' 


(38) 


It will be observed that for A < \/l3/3 there are two real roots. If 
A = 1, then f 1 / 3 = 1, 13/9. On the other hand, if A > \/l3/3, the roots 
are imaginary, so there are no intersections. If A = \/l3/3 there is a 
double root at t 1 / 3 = 13/11. In this case the parabolic trajectory just 
grazes the metric interface. 

In figure 1, I illustrate the two different trajectories when ro = 
\/l3/3 and <po = 1.0. In the case where ro is just any arbitrary amount 
smaller, the expanding interface overtakes the test particle following its 
Friedmann-Lemaitre parabolic trajectory and it must then follow the 
static Schwarzschild equations of motion. The Schwarzschild metric 
takes over at t = (13/11) 3 as explained above, and after that the 
trajectory is an ellipse with semimajor axis 7.6287 and the semiminor 
axis 3.986. These imply an eccentricity of 0.690 68 and a semilatus 
rectum parameter of p = 2.086. On the other hand, if ro is any arbi¬ 
trary amount larger, it escapes the moving interface and continues to 
follow the parabolic trajectory. It is evident from figure 1 that future 
trajectories are, in some cases, discontinuous functions of the initial 
conditions for the “Swiss cheese model.” 
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Put another way, the Schwarzschild metric permits closed orbits and 
the Friedmann-Lemaitre metric does not. In terms of the latter coordi¬ 
nates, one can choose initial conditions so that the parabolic trajectory 
just grazes the metric interface (fixed radial coordinate in this metric) 
and the speed is low enough that for infinitesimally different initial 
conditions the trajectory crosses the interface and is caught in a bound 
state, or alternatively misses the interface and proceeds on its parabolic 

trajectory. All these effects take place in supposedly empty space of the 
_2/3 

order of 25 h 50 million AU from a mass concentration of size Mq, and 
are quite counter to one’s physical intuition that such discontinuities 
should not occur there, except perhaps if the central void was created 
by an explosion leaving a remnant star behind. 

A third answer was provided by Bona and Stela[15]. It is a different 
type of “Swiss Cheese Model.” They insert in a flat space, Friedmann- 
Lemaitre-Robertson-Walker background a spherical patch in which, 
instead of the Schwarzschild metric, the Liebovitz[16] solution is in¬ 
serted. In the classification of Krasihski[10], this solution is in the “/T ^ 
0” subfamily of the Szekeres-Szafron family of solutions. Briefly put, 
this solution results from setting T j 1 equal to its Friedmann-Lemaitre 
value, i.e., by (9), 

( . \ 2 2 

-J -A = 2^ 1 +(—) -A (39) 
c I ac l \acJ 


Since the right hand side of (39) is independent of p, the spatial curva¬ 
ture is isotropic by (10), and is homogeneous. 

As (39) is a generalized Riccati equation for /3, it can be integrated. 
The general solution is, 


rr rlt 

e ^P’ T '> = a (r) [Ci(p) + C 2 (p)w(t )] 2/3 , where w(t)= -t-to- 

Jto «r) 

(40) 

For this solution, we can compute that, 


8 ttT 4 4 = 3 


ci 

- + 

a 


2 C 2 (p)w 

3[Ci(p) + C 2 {p)w\ 


(a 2C' 2 {p)w 

l« + 3[C[(p) + C' 2 (p)w\ 


(41) 


The dynamics of this model reveal that its predictions are at con¬ 
siderable variance with the observations of planetary astronomy. The 
slow-speed limit of the equations of motion follows directly from (29) 
and is, when we remember that Cq = 0(c -1 ) in the slow speed limit, 

a 2 C 2 {p) 2 _ 2aC 2 (p)r 

a 9a 3 r 2 3a(ra) 3 / 2 r 3 ’ 
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where r(p,r ) = e^^ p,T \ In order to agree with Newton’s result, we 
must choose C' 2 (p) = —§-\/2 GM. The choice of C!\ (p) does not affect 
the dynamics, so we pick C\ = ^^/2GM / c 2 p. This choice suffices to 
cause X 4 —► 3(a/a ) 2 as p — > 00 . (Note that different choices of C\ and 
C 2 are required for the zero mass case.) 

This result is substantially unique for a model which has isotropic 
spatial curvature, a non-zero central mass concentration, and is asymp¬ 
totic, at large distances from that mass, to the Friedmann-Lemaitre, 
expanding universe model. However, this model leads, by (42), to the 
results, 


a a [2GM 

~ r ^ 2 \— 

a a z V ar 


GM c 2 C 2 


(43) 


The problem to be noticed is the factor of a 3 in the denominator of 
the Newton gravitational attraction term. The import of this factor is 
that the radius of the earth’s orbit will be proportional to a 3 as the 
gravitational term and the angular momentum term dominate the other 
terms. By Kepler’s law that the ratio of the squares of the periods is 
proportional to the cube of the diameters, the year will be proportional 
to a 9 / 2 This means that the year should be increasing at the rate of 
about 0. 73/150 seconds per century. Here /150 is the Hubble constant 
in units of 50 kilometers per second per Megaparsec. Since the value 
currently accepted by the International Astronomical Union is +0.0095 
seconds per Julian century [17], the predictions of this model are in very 
serious disagreement with observations. The above results are in terms 
of the time variable r. Properly done, we need the results in terms of 
time measured on the Earth for the Liebovitz line element. The re¬ 
expression of this result in terms of “Earth time” involves a correction 
of the order of HqGMq /(c 2 R®) where M 0 is the mass of the Sun and 
r 0 is the radius of the orbit of the Earth. This quantity is of the order 
of 10 -8 smaller than Hq which is of the order of the predicted effect. 
Hence, this latter correction is negligible in this case. There is, of course, 
also a time independent correction to the relation between earth time 
and r of the order of GM & /(c 2 R 0 ) but it is not relevant to the effect 
discussed here. 

The conclusion of this section is that in one way or another, all the 
models reviewed that join the Schwarzschild metric at small distances 
and the expanding universe metric at large distances have problems of 
one sort or another. Panek [18] has stress the need to go beyond the 
“Swiss Cheese Model” and has analyzed numerically several different 
density profiles. In the next section, we will analyze the effect of a mass 
concentration inserted in a universe with a homogeneous mass-energy 
density. 


effects2r.tex; 1/02/2008; 20:58; p.12 



ACCELERATING EXPANSION OF THE UNIVERSE 


13 


4. HOMOGENEOUS, MASS-ENERGY DENSITY 

SOLUTIONS 


In line with the idea that the vacuum has a homogeneous self-energy, we 
consider homogeneous mass-energy density solutions to Einstein’s field 
equations. Onr case of interest is the embedding of a mass concentration 
in expanding, curved space. We impose the condition, 


8vrT| = 3 


i] + 2^- + e ~^ 

c / c 2 fi' 


1 ~f\p)~ 


2 f(p)f(p) 
& 


-A 


= 3 


d_y 3 
ca) (aRo) 2 


(44) 


which is the Friedmann-Lemaitre value in curved space. Since the right 
hand side is a function of r alone, we may integrate this equation, as 
at (15), with respect to p to give, 

P 2 + c 2 e _2/3 [1 - f 2 (p)] = (£) + + C 3 (T)e- 3 ^. (45) 


At short distances, we wish to match the Lemaitre-Schwarzschild metric 
(6). We find, upon substituting in (45), for 0 that 


2GMe~ 3f) + c 2 e~ 20 [l - f 2 {pj 



+ C 3 (r)e-^. (46) 


Thus, we shall choose C 3 (r) = 2GM independent of r. For later con¬ 
venience we rewrite (45) as 

P 2 = (£) + ^2 + 2 GMe-W - c 2 e~ 2/3 [l - f 2 (p)} . (47) 


If we are given /3(po, r o), then this non-linear, first-order differential 
equation may be integrated with respect to r to give P(po,r) V r. 

In the special case M = 0 we have the Friedmann-Lemaitre solution, 
and the special case, a = a = 0, Rq = oo gives us agreement with the 
Lemaitre form of the static Schwarzschild metric (6). 

Next we compute from (9) and the differentiation of (47). We 
obtain, 


8 tt c 2 Tl = 3 



3c 2 

(aR 0 ) 2 


(48) 
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\ 3 dc 2 

) « 3 ^o_ 

which reduces correctly in the two special cases “flat space” with d = 
a = 0, and M = 0. Since Tj 1 is not independent of p, the spatial curva¬ 
ture is neither homogeneous nor isotropic. It varies, as the distance from 
the mass concentration increases, from 3(a/a) 2 + 3c 2 /(ai?o) 2 according 
to (48), to the Friedmann-Lemaitre limit 2(a/a) + (d/a) 2 -|-c 2 /(ai?o) 2 as 
the distance tends to infinity. The McVittie solution mentioned above 
has isotropic, but inhomogeneous, spatial curvature. The difference 
with this solution lies in the choice here that the universal expansion 
factor is a function of the proper time and his choice that it is a function 
of his coordinate time, which is not an invariant quantity. 

So far, we have not had to discuss j3' or /?', as $ given by (47) 
has sufficed for all our computations. The general equations do not 
determine the behavior of [3 as a function of p. 

We observe that in the case of the expanding, curved-space of the 
Friedmann-Lemaitre model universe, that 

- cV 2 * 5 [l - f(p)] = 0 (49) 

Since we wish our solution to tend asymptotically to this limit, we 
choose to require this result in the long distance limit. As a result of 
this consideration, we choose as initial conditions for (47) 



ex P P(p,T 0 ) 

f(P ) 



(50) 

(51) 


where, as the dependence on Rq is actually a dependence on i?g, there 
is no problem with the case where < 0. The continuation is just 
.Tsin(p/a;) e-*• |x| sinh(p/|a;|). Using these initial conditions and a(r o) = 
1 we find that (49) is exactly satisfied initially. In addition the ratio 
between the surface area of a sphere and the square of the radius is in 
accord with that for the Friedmann-Lemaitre solution [5] for expanding, 
curved space in the open, closed, and “flat” cases. This choice, together 
with (47), defines (3 V p, r. 
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The dynamics are given by (29). The slow speed limit of the equation 
of motion is, 


f 



GM c 2 C 2 



2 GM 

3 


gg J 

^ [i - m 


(52) 


The curvature of space and the expansion of the universe contribute 
corrections of second order which are currently undetectable in the 
planetary motions of our solar system. 

The solutions of this section have assumed that: (i) There is spherical 
symmetry about the mass-energy concentration, (ii) Tf = 1 4 = 0. (iii) 
is homogeneous as specified by (44), but may be time dependent, 
(iv) The metric matches in the short distance limit, the form of the 
Lemaitre form of the static Schwarzschild metric corresponding to a 
mass-energy concentration of strength M. That is to say, the region in 
which the gravitational force is much greater than that of the expansion 
of the universe, (v) The metric matches the long distance form of the 
Friedmann-Lemaitre metric. Under these assumptions, the solutions 
of this section is substantially unique, subject to the usual freedom 
involving changes of variables. This result is a parallel to Birkhoff’s 
theorem[19], but involves a different class of models than the class he 
considered. 


5. PROPERTIES OF THE HOMOGENEOUS, 
MASS-ENERGY DENSITY SOLUTIONS 


It is convenient to make the change of variables, 

A T) = o(r)e ffW . 

Then (47) becomes, 


(53) 


• 2 „ u 0 


+ 2 GMe~ 30 - c 2 e~ 2f) [l - f 2 (p)} , 


which may be solved by means of the quadratic formula to give 


a _i_ 

£7 = — ± 

a 


+ 


(■aRof 


+ 


2 GM 




c 2 p 2 
R 2 e 2 0 


11/2 


(54) 


(55) 
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where I have chosen a simple reparameterization of (45-51), to wit, 


p = R 0 sin (56) 

so that (51) becomes, 

f(p) = \/ 1 - p 2 / r o ( 57 ) 

I choose as initial conditions for (55), 

a(r 0 ) = 1.0, exp [a(p,T 0 )\ = p (58) 


We may now check the special case, M = 0 and Rq = oo. As we 
expect = a(r)p, it follows that a = 0 for this case. That result in 
turn implies the minus sign in (55) as a is known to be positive. Next 
we check the case where e^ is very small. The dominant terms in (55) 
are, 


& = ± 


2GM 


(a(r)e' 


<r\3 ' 


(59) 


Integrating this equation with respect to r yields, 


o / 2 C A/f C T 

exp [§ <r(p,T)\ - exp [§ cr(p,T 0 )\ = c J a~ 3/2 (t)dt (60) 

When we note that the integral in (60) is approximately t — To, and by 
use of (58), we get, 


exp [^(j(p,T) 



2GM . 

— 5 —c(r - r 0 ) 


(61) 


In order to match the Schwarzschild-Lemaitre metric ( 6 ) for small 
we must choose the minus sign in (61) and hence also in (55). From the 
examination of these special cases, we conclude in general (assuming 
that a > 0) that (55) must always be taken as, 


U = —U-- 
2 a 2 



2 GM 
a 3 (r) 


+ 


c 2 VlP 

(a(r)Ro) 2 



(62) 


where U = exp(|cr). 

It is worth pointing out that when a = 1.0 that (62) subject to the 
initial conditions (57-58) gives the static, curved-space solution to the 
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Einstein field equations for a mass concentration in a vacuum. In this 
special case, P(p,r) = a(p,r) by (53). Eq. (62) becomes in this case, 


U = — 


3 

2 


2 GM + 



(63) 


where here U = exp(|/3). The integration of this equation follows 
directly using the aforementioned initial conditions. a(p, r) follows by 
( 8 ) and (57). By (44) and (48), we conclude that the T 4 = Xf = X 3 = 
3/(8ttRq) and 871 X 4 = Rq 2 , as expected in a curved-space vacuum. 

The time when an observer at rest with respect to the coordinate 
system, i.e. a freely falling observer, arrives at the mass concentration, 
i.e., when e a = 0 for a given p with initial conditions (58), is given by 


U(p,T zero ) - U{p,T 0 ) = - p 3/2 = 


3 [ Tzei ° , / [ 2 GM fa\ 2 C 2 ^IP 
2 J tq \ a 3 (r) + \a ) + (a(r)X 0 ) 2 



2 p 3 / 2 

Tzer ° “ T ° + 3^2 GM 



1/2 

(64) 


by (61) and (62) Thus the physically allowed range of time is 0 < r < 

Xero ip) ■ 

For numerical computations it is convenient to choose dimensionless 
quantities. I select, 


X = H 0 t, U 


Ro 




U, M 


GMHq A = ^_ Ro. 
c 3 ’ A H$a c P' 


(65) 


The initial conditions are ^4(X 0 ) = A(Tq)/A{T q ) = 1,11 (R, X 0 ) = T7 3 / 2 , 
and TWq ~ 0.799 x 10 ~ 23 /i 5 o. For reference, one astronomical unit is 
about 0.808 x 10 _ 15 /i 5 o in our dimensionless units and a Megaparsec is 
about 1.67 x 10~ 4 h^. Thus we may rewrite (62) as 


U = 3 4u- 3 - 

2A 2 


A V 2 M y/U 2 

A I + A 3 (T) + (A(T)R 0 ) 2 


1/2 


(68) 

I have integrated this equation by means of the Runge-Kutta method 
[ 20 ] adapted to double precision. 
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There are several special cases for the solution of equation ( 66 ). 
First, when is large, the solution is just = pa(r), both for the 
case of “flat” and curved space. This solution is just what one expects 
for the case of expanding Friedmann-Lemaitre expanding space. The 
quantity e a is given by (57) and ( 8 ) which is, of course, different for 
“flat” and curved space. I illustrate the large case, and the approach 
to it in fig. 2 for the case of flat space and a(r ) = ( t/to ) 2 / 3 . The initial 
conditions (58) insure that e 13 /\pa(ro)\ = 1 in every case. The central 
mass has been chosen to be that of the Sun. 

The next special case is when is small. In this case, in line with 
with the Lemaitre form of the Schwarzschild metric, it is appropriate 
to plot the results vs. (r — To)/p 3 A This is because this form of the 
(static) metric can be written as 


eP_ 

P 


3 /2 GM f c(r-r 0 )\ 

2 V c 2 V P 3/2 J 


(67) 


I illustrate this case and the the approach to this limit in fig. 3. The 
large dots indicate the expected [p — c(r — To )] 2 / 3 limiting curve with 
vertex at r zero . 

The divergence at r = 0 noted in fig. 2, leads to the question of the 
nature of this divergence. It turns out that it depends on the exponent 
7 defined in the limit as r —► 0 by a(r ) x IFr 7 . By considering the 
dominant terms in ( 66 ) for flat space, we obtain the results, 


e ' 3 


c 

H~o 

c 

#0 

c 


H 0 


2M 


4( 7 -i ) 2 - 7 2 


1/3 


^y /3 

T o) 


7 r 


Cl{p) ~ 


T Y 


3 M 


7(2 - 3 7 ) 


r 2-3 7 


1 1/3 


--M — In — 


r o J 


TO 


1/3 


7 > |» ( 68 ) 

0 < 7 < |(69) 
7 = | (70) 


The behavior for values of 7 not illustrated is qualitatively similar, 
but differs, of course, in detail. 

The effects of the curvature of space on (3 is not very significant 
either for small values of e^ or for large values, as in these cases the 
limiting values discussed above are obtained. The effects however are 
of some significance for intermediate values of e^. I have illustrated an 
example in fig. 4. 
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6. DOMAINS OF ATTRACTION 


In this section we will focus our attention on the slow-speed (52), flat- 
space[l, 2] limit of (29) for radial motion only. That is to say, Cq = 0. 
Hence, using the dimensionless variables of (65), we obtain, 



where 91 = /c. First consider the case where the rate of expansion 

of the universe is decreasing. Then A < 0. If we divide (71) by (A/A) 2 91 
and let x = MA 2 /(91 3 A 2 ), then it becomes 


fA\ 2 91 _ 1 _ ^ _ 1 ^ AA 

~~ X \/l + 2x + A 2 A l + 2x 


(72) 


Note that the A term in (72) is uniformly negative over the allowed 
range of x, (0 < x < oo). Hence we may check all the cases A < 0 
by setting it to zero. This being done, we may compute in a straight 
forward manner that for very small x, the right-hand side is just —1.5.x 2 . 
Differentiation with respect to x shows that the derivative is uniformly 
negative. Thus we conclude that if the rate of expansion of the universe 
is not increasing, then according to the metric of found in section IV, 
there is just one domain of attraction and it is centered on the mass 
concentration as (72) shows that if A < 0, then 91 < 0. The analogous 
computations for curved space proceed in a similar manner, but are a 
bit more complicated. 

Next let us consider the case where, as reported by Reiss et al. [3] and 
Perlmutter et al. [4], the rate of expansion of the universe is increasing. 
As for this case, A > 0, it follows directly from (71) that in the limit 
as 91 —» oo, 91 > 0. Since it is also the case that when 91 —> 0, 91 < 0, 
it follows that there exists an 91 o for which 91 = 0. Therefore, we now 
have two domains of attraction. One is 0 < 91 < 91In this domain, 
all test particles at rest with respect to the coordinate system appear 
to be accelerating towards the origin. The other domain is 91 d < 91. 
In this domain all test particles at rest with respect to the coordinate 
system appear to be accelerating towards infinity. One expects that 
over time, the neighborhood of the domain boundary will be swept 
clean of particles. 
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Note is taken that the “Swiss cheese model” also has 2 domains of 
attraction. This fact is independent of the sign of a. Also the location 
of the boundary described by IK increases with time. 

At this point it is convenient to define the standard acceleration/de¬ 
celeration parameter qo = — AA/(A 2 ). Then, setting Tt = 0, the equa¬ 
tion for the domain boundary becomes, by (72), 


0 = 1 — x — 


Qo + 1 
n/1 + 2x 


(73) 


The value of qo is given[21] by the equation, 

Qo = ~ Ha (74) 


where, combining the results of [1, 2, 3, 4] we deduce the values Qm = 
8itGpo/(3Hq) & 0.19 and 12a = Ac 2 /(377q) 0.81. These values yield 

q o ~ —0.715. Note that for “flat” space, —1 < qo, as 0 < Qm- Some 
simple manipulations convert (73) into a cubic equation in x. 

2 x\) — SxJj + 1 = (1 + qo) 2 (75) 


This equation is simply solved by standard methods. Its solutions, of 
course, also include those for the other sign of the square root in (73). 
Evidently, the solutions of (73) all lie in the range 0 < x < 1. If qo > 0, 
there are no positive real roots of (73) and we take xd{qo) = 0. If 
qo < 0, then there is one positive real solution xd{qo) of (73). For the 
above mentioned case, £r>(—0.715) « 0.825. 

We now have the domain-size-mass relation, 

(76) 


r D = 


GMa 2 
\x D (q 0 )d 2 \ 


Here ru is the radius of the domain and M is the interior mass. We 
appeal to the uniqueness (in our class of models) reported at the end 
of the fourth section to allow us to treat all the mass as concentrated 
in the center of the domain. See also the result (16) in this regard. 
The domain-size-mass relation is an important conclusion which is of 
particular interest if the rate of expansion of the universe is accelerating, 
as has recently been reported[3, 4], 

It is to be noted that many mass concentrations are rotating and so 
are only axisymmetric, rather than spherically symmetric. This means 
that near the mass concentration the metric has some of the charac¬ 
teristics of the Kerr metric. In terms of Cartesian coordinates, the dif¬ 
ferences in the far-held between the Kerr metric and the Schwarzschild 
metric [8] are x 2 G(r x I ■ ei)/(rc) 3 , where the e) are unit vectors 
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in the x, y, z directions and I is the angular momentum of the central 
mass concentration. In regions where this correction is relatively small 
compared to the gravitation term 2 GM/(c 2 r), i.e., \r x I\/(Mcr 2 ) <C 1, 
we may safely neglect this effect. On the other hand, where the correc¬ 
tion is significant, the precise details of this effect will require further 
study. In addition, our lack of a clear knowledge of the distribution of 
dark matter also creates uncertainty concerning the details the above 
obtained results. 

We now consider a few applications of the domain-size-mass rela¬ 
tion. If we use the solar mass, we get the radius of the domain boundary 
—2/3 

to be about 417 h 50 light years. However, in this case the density of 
stars is sufficiently high so as to perturb significantly this result. The 
Oort cloud only extends out to about 0.8 light years. We remark that 
the domain-size-mass relation only provides an upper limit on the size 
of a structure that a given mass can be expected to hold together. That 
there are smaller structures can presumably be explained on the basis 
of other circumstances. 

Next consider a globular cluster. They are observed to have about 
10,000 to 1,000,000 stars and to have diameters of several 10’s of 

light years to about 200 light years. The computed domain boundary 

—2/3 —2/3 

diameters for systems of this size are about 9000 h 50 ' - 42000/r 50 

light years. This result is of the substantially larger than that which is 
observed. It may be that they are reduced in size when they orbit near 
the galactic center. 

Let us now consider our galaxy. It is estimated to contain something 
like 200 billion stars and to have a mass of 700 billion to 1 trillion solar 
masses. This mass implies a domain radius of 3.7/i 50 ^ 3 to 4.2 h 5 ^ 3 
million light years, which way off from the observationally estimated 
diameter of about 100 thousand light years. 

When we consider the Local Group and the Virgo Cluster, the 
problem becomes more interesting. Specifically, it has been estimated 
that the roughly 30 galaxies in the local group have a mass of about 

3 x 10 12 Mq, not counting intergalactic dark matter. This mass corre- 

—2/3 

sponds to about a radius of 6 /i 50 ' million light years. The observed 
value is roughly 4 million light years. In the Virgo Cluster of around 
2000 galaxies, the mass is thought to be of the order of 10 15 Mq. This 

mass includes a significant amount of dark matter and gives a radius 
—2/3 

of about 42 /i 50 w million light years, compared to the observation that 
the Milky Way galaxy lies about 65 million light years from center of 
the cluster. Considering the errors and the fact that the Virgo Cluster 
appears not to be in equilibrium yet, this agreement is not too bad. 
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In conclusion, it appears that for larger scale structures composed 
of galaxies and inter-galactic space, the observed increase in the rate of 
expansion may be an important feature in determining the size of self¬ 
bound gravitating systems. For smaller structures like galaxies, globular 
clusters etc. other mechanisms are presumably dominant. 
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Figure 1. “Swiss cheese model” trajectories which begin at xo = \/l3/3 + e (the 
Friedman-Lemaitre case), and begin at xo = \/l3/3 — e (the Schwarzschild case). 
Here, e > 0 may be chosen as small as one pleases. The initial part of both 
trajectories is a parabola generated by the Friedmann-Lemaitre metric. At time 
t = (13/11) 3 , marked in the figure, the two trajectories separate. 
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Figure 2. The behavior of e /3< - p, ' r ' ) for larger values of p. The solid line is the limiting 
value obtained for large values of p. It is e 13 = a(r)p. The other lines show how the 
limit is approached as p increases. 
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Figure 3. The behavior of e l3( - p,T ' > for smaller values of p. The solid line is the limiting 
value obtained for smaller values of p. It is e^ 3 oc a(r)p[l — c(r — T zero )/p 3 / 2 ] 2 / 3 . The 
large dots show this formula. The other lines show how the limit is approached as p 
decreases. 
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r - r (units of the Hubble time) 


Figure 4- The variation of e /3 ^ p ’ r ' ) with the curvature of space. Displayed are the 
ratios of e 13 for curved space to that for “flat” space with the same values of p, r. 
The sold curve is for Rq = +1 and the dashed curve is for Rq = — 1. 
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